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ABBREVIATIONS 


The  following  abbreviations  are  used  in  this  document: 


ADINA 


-  Automatic  Dynamic  Incremental  Nonlinear  Analysis  Code 


CAMPACT  -  Container  for  Armaments  Protection  and  Transport 
DOF  -  degree(s)  of  freedom 

FE  -  finite  element 


'inite  element  method 


-  H&R  Technical  Associates 


115  mm  chemical  agent  rocket  weapon 


Southwest  Research  Institute 


-  two  dimensional 


-  three  dimensional 


1.0  BACKGROUND  AND  SUMMARY 


H&R  Technical  Associates  r.as  been  contracted  to  assess  the  risks 
associated  with  the  transport  of  M55  chemical  agent  munitions.  In  support  of 
these  studies,  SwRI  has  analyzed  for  H&R  various  sub-problems  in  the  areas  of 
thermal  cook-off,  structural  damage,  and  penetration  of  these  munitions  and 
their  shipping  containers.  Separate  reports  have  been  issued  by  SwRI  stating 
the  results  of  these  analyses. 

This  particular  document  addresses  the  impact  analyses  performed  on  the 
rocket  packing  assembly,  hereinafter  referred  to  as  the  "pallet  assemoly." 

The  concern  motivating  these  analyses  can  be  stated  as  follows.  What  is  the 
highest  distance  from  which  the  pallet  can  oe  dropped  onto  an  unyielding 
surface  without  inducing  leakage  in  in  the  agent  cannister?  In  response  to 
this  concern,  a  pallet  of  M55  rockets,  illustrated  in  Figure  i,  was  analyzed 
to  determine  if  it  could  be  dropped  onto  a  unyielding  surface  from  40  feet  on 
each  of  its  faces  (end,  side,  and  bottom  faces)  and  not  induce  failure  of  the 
agent  cannister. 


It  is  anticipated  tha 
in  a  protective  containers 
resemble  ordinary  20-foot 
constructed  for  transporti 
construction  features  incl 
and  Kevlar  sidewalls,  and 
truss,  and  inner  lining  pr 
thermal/fire  protection  ov 
illustrates  this  container 


t  pallet  assemblies  of  M55  munitions  will  be  shipped 
named  "CAMPACT."  On  the  exterior,  CAMPACTs 
shipping  containers.  They  are  especially 
ng  hazardous  materials,  however.  CAMPACT 
ude  a  thick  foam  inner  lining,  aluminum  honeycomb 
an  interior  truss  for  added  strength.  The  sidewall, 
ovide  a  substantial  increase  in  rigidity  and 
er  conventional  container  construction.  Figure  2 


If  a  loaded  CAMPACT  container  is  dropped,  the  foam  lining  of  the  CAMPACT 
might  have  an  ameliorating  effect  on  the  agent  cannister  response.  On  the 
other  hand,  the  additional  mass  of  pallet  assemblies  above  the  agent  cannister 
of  interest  might  significantly  increase  the  cannister  response.  To  establish 
in  a  qualitative  manner  the  effects  of  the  foam  lining  and  pallet  interaction, 
additional  calculations  were  also  performed  on  a  simpie  pallet  model  modified 
to  reflect  the  cushioning  of  the  foam  inner  lining  and  the  added  inertia  of 
other  pallets.  (For  detailed  calculations  of  the  structural  response  of  the 
CAMPACT  to  impact  see  [ 1 ] ) . 

It  is  important  to  briefly  mention  at  this  point  some  conditions  that 
were  assumed  to  make  the  impact  analyses  tractable  and  the  scope  of  the  work 
performed.  Firstly,  it  was  assumed  that  the  agent  cannister  was  in  a  new 
condition,  that  is,  aging  effects  such  as  load  history,  corrosion,  creep  and 
pallet  degradation  were  not  in  any  way  assessed  or  accounted  for.  Secondly, 
it  was  assumed  that  agent  cannisters  were  not  leaking  at  the  time  of 
transportation.  Thirdly,  in  regards  to  the  scope  of  this  work,  these  analyses 
quantitatively  address  only  the  impact  of  the  pallet  on  a  unyielding  surface, 
separate  from  any  shipping  container  which  might  be  used  to  transport 
pallets.  Interaction  effects  nave  been  addressed,  in  a  qualitative  way,  for 
the  longitudinal  impact  case,  nowever. 

Results  of  the  pallet  impact  analyses  can  be  briefly  summarized.  In  the 
case  of  bottom  or  vertical  impacts  onto  unyielding  surfaces  from  40  foot 
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heights,  it  was  concluded  that 
occur.  For  end  or  longitudinal 
analyses  indicated  that  strains 
criteria  and  leakage  of  the  age 
Longitudinal  impacts  onto  unyie 
failure  of  the  agent  cannister, 
of  the  pallet  onto  unyielding  s 
failures  would  probably  occur, 
pallets  and  the  CAMPACT  indicat 
ameliorating  effect  on  tne  pal- 
foot  heights. 


failure  of  the  agent  cannister  would  not 
unyielding  surface  impacts  from  40  feet,  the 
in  the  agent  cannister  exceeded  the  failure 
nt  cannister  probably  would  occur. 

Idir.g  surfaces  from  30  feet  did  not  cause 
however.  Analyses  of  sice  or  lateral  impacts 
ur faces  from  40  feet  showed  that  cannister 
Calculations  of  the  interactions  between 
ed  that  the  CAMPACT  could  have  a  significant 
et  response  during  longitudinal  impacts  from  40 


This  report  has  been  organized  to  assist  the  technical  reader  in 
understanding  the  results  achieved,  and  also  how  they  were  derived.  It 
consists  of  two  major  parts.  The  first,  or  methodology  section,  describes  the 
problems  solved,  the  details  of  the  finite  element  models  used  to  generate 
strain-time  histories  in  the  agent  cannister,  and  general  features  of  the 
nonlinear  computational  scheme.  The  second  part  presents  the  results, 
primarily  in  the  form  of  strain-time  histories  for  nodes  at  several  positions 
along  the  agent  cannister.  These  plots  are  discussed  in  light  of  the  failure 
criteria  adopted  and  conclusions  are  then  drawn  regarding  leakage  of  agent. 
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ANALYSIS  METHODOLOGY 
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It  does  rot  require  tiucr.  stretching  of  the  imagination  to  suspect  that 
cppi-g  pallets  frem  -0-foot  heights  onto  unyielding  surfaces  could  cause 
e.di-.g  cr  failure  of  pallet  components.  Responses  of  this  type  nave  the 
r.erai  mass  if  i  cat  ion  of  nonlinear  response,  ’//hat  is  meant  by  nonlinear 
scorse  is  tr.at  computed  displacements  are  not  linear  functions  of  the 
p_.ee  .cads.  This  is  a  mathematically  correct,  but  physically  an 
appeal.  1  mg  explanation. 


Several  physica.  phenomena  can  oe  identified  as  causing  nonlinear 
spon.se:  material  yielding,  changes  in  boundary  conditions,  and  large 
sc .acements .  Material  yielding  resuits  in  a  nonlinear  relation  cetween 
ter. a.  stress  and  strains:  large  displacements  can  cause  yielding  also,  but 
ey  crar.ge  the  directions  of  the  applied  loads  on  the  structure  as  well., 
arges  in  ro-r.cary  conditions  generally  occur  when  the  structure  makes 
ntac:  -ith  otr.er  structures,  usually  after  some  initial  displacement.  This 
tier  resu.ts  in  a  nonlinear  response  oecause  there  is  a  "step  change"  in 


A-alvses  of  the  pal.et  assembly  involved  all  these  types  of 
"  .  .rear  :  t ;es.  It  was  nearly  certain  from  the  outset  of  this  project  that 
e  att.mc  .sacs  ar.d  low  yield  points  cf  some  pallet  components  would  result 
.re.ast.:  materia,  cenavior.  Although  not  so  obvious,  it  also  became  clear 
at  .arge  displacements  would  be  induced  as  well;  seme  component 
sc. acements  exceeded  one  inch.  Finally,  the  structure  of  the  pallet  is  suen 
it  - : 1 1  large  displacements,  some  components  could  make  contact  with  the 
tact  s-rface  and  therefore,  changes  in  boundary  conditions  could  also  occur. 


?r  program  oreposed  and  used  for  the  pallet  impact  proolems  was 
i.  four  levels  of  computational  complexity  exist:  linear 
illy  nonlinear,  total  Lagrangian,  and  updated  Lagrar.gian 
computational  .evel  is  more  complex  because  fewer  assumptions 
isral  response  are  made,  e.g.,  small  strains. 1  inear  materials, 
arms.  ar.d  constant  boundary  conditions.  It  appeared  at  the 
:'ew  assumptions  like  these  could  be  made  in  the  case  of  pallet 
5  resu.ts  have  indeed  validated  this  approach.  Thus,  for  the 
'oolems,  an  updated  Lagrangian;  the  most  computationally 
Lque,  was  employed  to  construct  solutions.  That  is,  materials 
oehave  nonilrearly,  displacements  were  allowed  to  become  large 
loading  directions  and  boundary  conditions  were  accounted  for. 


-as  ceer  minted  out  that  one  of  the  features  of  these  impact  problems 
:ert  am  assumptions  male  in  more  common  engineering  problems  could  not 
-.trc.gr  it  was  not  possible  to  reduce  the  computational  effort 
i  to  ctr.struct  solutions,  it  was  possible  to  reduce  the  problem  size, 
e.ement  motels  ~ere  kept  to  a  reasonable  size  (less  than  250 
tec  is  equations)  by  making  plausible  assumptions  about  pallet  behavior 
i  tartic-.ar  impact  orientation.  These  presumptions  about  pallet 
response  set  tr.e  general  features  of  the  FE  models  that  were 


1 . i  leve i oted 
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First,  it  was  presumed  that  in  each  of  the  various  impact  scenarios,  only 
certain  elements  of  the  pallet  assembly  would  te  significantly  loaded.  Thus, 
the  rest  of  the  pallet  need  not  be  explicitly  modeled  so  long  as  its  inertia 
was  accounted  for.  A  second  presumption  made  *as  that  the  impact  response 
consisted  of  rocket  bending  and  cross  section  crushing  loads  (at  the  rocket 
supports)  and  that  these  phenomena  could  be  treated  separately.  Thus,  ber.aing 
and  crushing  strains  were  computed  using  different  FE  models  ar.d  results  were 
superposed . 


In  consequence  of  these  initial  assumptions  about  the  pallet  resoon.se 
behavior,  five  finite  element  models  were  created  to  analyze  bottom 
(vertical),  side  (lateral),  and  end  (longitudinal)  face  impacts.  Vertical  and 
lateral  impact  analyses  utilized  two  models  each:  a  beam  ar.d  truss  element 
model  of  one-quarter  of  the  pallet  assembly  to  obtain  oer.dir.g  strains,  and  a 
two-dimensional  model  of  the  rocket  support  region  to  ootain  strains  in  the 
agent  cannister  due  to  crushing  or  squeezing.  The  end  face  or  longitudinal 
impact  scenario  does  not  involve  any  crushing  or  squeezing  of  the  agent 
cannister  between  lateral  suooorts  so  on^v  one  model  was  constructed. 


Assumptions  about  pallet  response  oehavior  also  guided  the  selection  of 
the  nonlinear  material  models  used  for  pallet  components.  Where  strains  were 
expected  to  remain  in  the  elastic  range  linear,  elastic  models  were  used. 

Where  strains  might  exceed  yield  values,  elastic-plastic  models  were  selected 
appropriate  for  tne  material.  Nonlinear  material  models  used  in  these 
analyses  were  isotropic  elastic-plastic,  elastic-perfectly  plastic,  nonlinear 
elastic,  or  orthotropic  elastic  [3].  For  nonlinear  elastic-plastic  models,  a 
material  elastic  modulus,  yield  strength,  and  tangen:  modulus  are  specified; 
isotropic  strain  hardening  is  assumed  in  this  material  model.  Elastic- 
perfectly  plastic  models  appear  to  be  nearly  mathematically  identical  to 
elastic-plastic  models  but  the  tangent  modulus  is  input  as  zero.  These  models 
constituted  the  principal  material  representations  in  the  3-D  analysis. 
Nonlinear  elastic  material  models  were  invoked  to  represent  rocket  tube 
crushing  phenomena  and  the  CAMPACT  foam;  orthotropic  elastic  mater. ai 
were  used  to  represent  the  wood  in  the  rocket, tube  2-0  crushing  analyses. 


In  the  version  of  ADINA  available  at  SwRI,  there  is  currently  only  or.-: 
orthotropic  material  model  active:  a  linear  elastic  material  model.  For  tne 
2-D  models  described  below,  such  a  material  model  was  considered  acceptable. 
The  2-D  models  were  used  to  obtain  crushing  or  squeezing  stress-strain 
characteristics  of  M55  assemciies  resting  on  lateral  supports.  Only  localized 
plasticity  of  the  wood  supports  was  expected  in  these  analyses,  hence,  a 
purely  elastic  representation  of  the  support  was  considered  to  give  reasonable 
results.  In  the  3-D  mode.s,  on  the  other  hand,  large  amounts  of  plasticity 
were  expected  in  the  bending  of  the  wood  supports.  Thus,  a  linear-elastic 
model  was  not  sufficiently  accurate.  In  this  case,  isotropic  material  models 
were  used.  Material  properties  selected  for  these  "isotropic"  wood  materials 
were  the  orthotropic  parameters  corresponding  to  the  most  heavily  loaded 
orthotropic  axis.  For  example,  in  the  vertical  impact  model,  where  the  wooden 
lateral  supports  undergo  bending  about  axes  normal  to  the  grain  direction, 
grain  direction  mechanical  properties  were  input  as  the  "isotropic"  material 
model.,  as  the  principal  bending  stresses  would  occur  in  the  grain  direction. 


Because  the  version  of  ADINA  available  at  SwRI  does  not  contain  gap 
elements  -  elements  specifically  designed  to  address  contact  phenomena  -  it 
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was  difficult  to  adequately  address  tr.ese  conditions  in  the  lateral  pallet 
impact  proolem.  A  nonlinear  elastic  e.emer.t  with  increasing  stiffness  was 
created  for  the  vertical  pa. let  proolem;  and  used  with  some  success  to 
represent  contact  between  lateral  supports  ar.a  the  ground,  fetalis  are  given 
in  Section  2.4. 

Seme  assumptions  were  also  made  concerning  failure.  Failure  oy  excessive 
straining  was  assumed  to  occur  when  tr.e  total  strain  at  any  point  the  agent 
cannister  exceeded  a  predetermined  value.  The  Structural  Alloys  Hancboox 
gives  tr.e  ultimate  strain,  of  out ’-16.  tr.e  cannister  material,  as  3"  '  4  ] . 

While  a  strain  of  3?  at  failure  migr.t  ce  reached  in  a  test  coupon .  it  was 
believed  that  was  a  more  realistic  failure  strain  for  tne  agent 
cannister.  This  reduced  ai.cwao.e  strain  reflects  uncertainties  accut  the 
actual  agent  cannister  strength  capacity  ar.d  the  probability  tr.at  seme  stress 
concentrations  exist  for  the  cannister  which  are  not  accounted  for  ir.  the  -E 
model . 

Buckling  failure  modes  were  considered  for  er.d  face,  or  longitudinal 
impacts.  Failure  oy  buckling  was  determined  by  comparing  tr.e  critical  stress 
for  an  empty  cylindrical  shell  raving  cross-sectional  dimensions  equal,  to  tr.e 
agent  cannister)  with  the  maximum  ouexlirg  stress  caused  by  a  longitudinal 
impact.  Failure  by  buckling  was  presumed  if.  during  tr.e  impact,  loads 
exceeded  the  buckling  critical  value.  Mote  that  tr.is  approach  may  be 
conservative.  In  general,  tr.e  agent  cannister  is  considerably  more  than  half 
full  of  agent;  tne  fluid  ccu^d  sign  if  icantly  stiffen  the  agent  cannister  and 
thereby  increase  the  critical  buckling  stress. 

In  the  sections  below,  each  of  these  finite  element  models  are  discussed 
in  detail.  These  discussions  include  descriptions  of  the  nonlinear  material 
model,  active  degrees  of  freedom,  special  purpose  elements,  boundary 
conditions  assumed  and  the  applied  initial  conditions.  Models  are  discussed 
in  order  of  increasing  complexity:  longitudinal,  lateral  ar.d  vertical  impact 
analyses . 

2 . 2  Longitudinal  Impact 

2.2.’  Analyses  Performed 

End  face  or  longitudinal  impacts  of  pallet  assemblies  actually  involved 
analysis  of  only  the  M55  rocket  sub-assembly.  It  was  assumed  for  this  loading 
scenario  that  no  interaction  took  place  oetween  wooden  pallet  members  or 
fiberglass  launching  tubes  and  the  rockets.  (Calculations  in  Appendix  D  show 
that  this  assumption  appears  to  be  a  valid  one).  Therefore,  the  finite 
element  models  of  the  "pallet  assembly"  consisted  of  a  single  M?5  rocket, 
oriented  nose  downward,  with  boundary  conditions  simulating  a  rigid  impact 
surface.  Neither  the  compliance  of  the  plywood  end  caps  \see  Figure  1)  or  the 
launching  tube  end  plugs  were  included  ir.  the  longitudinal  impact  analyses. 
This  approach  is  conservative  although  it  was  judged  that  their  effect  on  the 
agent  cannister  response  wouid  ce  small. 

SwRI  numerically  calculated  the  nonlinear  dynamic  response  of  a  single 
F  rocket  impacting  unyielding  surfaces  from  heights  of  30  and  40  feet  and 
mated  the  effects  of  the  CAMPACT  structure  and  its  internal  pallet 
ai .  ^ment  on  the  safe/not  safe  drop  height. 
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Berore  computing  the  dynamic  impact  response  o:  the  M5t>.  the  finite 
element  model  was  checked.  Hand  oalouiat ions  were  made  of  the  rocket  weight 
and  the  first  longitudinal  natural  frequency.  The  estimated  rocket  weignt  of 
57  lbs.  agreed  very  well  with,  the  numerically  computed  weight  of  57.01  lbs. 

The  computed  first  longitudiral  frequency  was  210  Hz.  Closed  form  solutions, 
assuming  a  uniform  distribution  of  mass  and  stiffness,  gave  the  first 
longitudinal  natural  frequency  as  2"6.4  Hz.  These  closed  for-  solutions  are 
contained  in  Appendix  A. 

2.2.2  Model  Development 

To  analyze  impact  on  unyielding  surfaces,  the  investigators  ..set  a  oeam 
element  model  of  the  M55 .  This  model  consisted  of  28  beam  elements  having 
only  one  translational  degree  of  freedom  (in  the  rocket's  axial  direction)  at 
each  node.  Features  of  this  model  are  shown  in  Figure  3-  In  general,  only 
the  rocxet  casing,  anc  not  rocket  internals,  was  considered  in  establishing 
the  model's  axial  stiffness.  An  exception  occurred  ir.  the  warhead  region, 
where  effective  cross  sectional  procerties  were  computed  that  accounted  for 
the  additional  stiffness  o:  the  burster  casing.  Details  are  in  Appendix  A. 
Note  that  the  ourster's  internal  steel  sleeve  and  plastic  tube  w ere  .not 
considered  to  add  to  the  warhead's  axial  stiffness,  primarily  because  cf  their 
lack  of  longitudinal  end-fixity  15]. 

Since  tapered  team  elements  are  r.ot  available  in  ADINA,  tapered  sections 
of  the  rocxet  were  represented  cy  snort  beam  elements,  each  of  uniform  cross 
section,  out  increasing  in  area  in  the  direction  of  increasing  section 
diameter.  Each  element  in  a  tapered  region  had  outs:  le  diameters  equal  to  the 
average  outer  diameter  of  the  tapered  section  represented. 

To  achieve  a  correct  distribution  of  mass  along  the  rocket's  length, 
effective  densities  were  ccmouted  for  each  element.  These  effective  tensities 
accounted  for  rocket  casing  and  internal  weights  such  as  the  fuze  mecr.anism, 
the  chemical  agent,  and  the  solid  fuel.  The  exact  density  of  tr.e  chemical 
agent  could  r.ot  oe  determine:,  so  the  density  of  water  was  iset.  Because  the 
agent  density  was  unknown,  tr.ere  may  be  some  error  in  the  FE  model's  mass 
distribution.  Total  rocket  mass  used  in  the  analyses  was  correct,  however, 
and  equaled  57  ics. 

Each  M55  m  tr.e  pallet  assembly  is  ir.side  a  fiberglass  iauncr.ing  tube 
having  closed  ends.  These  t  .be  end  plugs  are  fabricated  from  aluminum  and  are 
counteroored  or  miilea  to  accept  the  rocket  tip  or  fir.s.  For  the  longitudinal 
impact  analysis,  the  plug  at  tr.e  nose  of  the  rocket  is  of  interest.  This 
fitting,  three  inches  thick,  is  counteroored  for  three-quarters  of  its 
thickness  so  as  to  accept  the  M55  fuze  assemoly.  As  this  puts  very  little 
material  between  tr.e  rocket  rose  and  the  impact  surface,  the  erd  plug 
stiffness  was  ignored  in  tr.e  calculation  of  rocket  impact  response.  Shear 
loads  between  the  cour.terbcre  in  tne  plug  ant  the  fuze  were  r.ot  computed. 

In  addition,  a  plywood  end  cao  at  each  end  of  the  pallet  prevents 
longitudinal  motion  of  the  rocket 'tube  assemolies  during  transport.  Erd  caps 
are  3.-  k-inch  thick  plywood  sheets  with  holes  located  at  each  launching  tuce 
end  plug.  End  cap  hole  diameters  are  smaller  than  er.d  plugs.  Again, 
compliance  of  the  plywood  end  cap  was  ignored,  resulting  in  a  conservative 
calculation  of  the  rocket  response  luring  longitudinal  impact. 
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Figure  3  also  snows  the  model  that  .,as  developed  to  estimate  the  effect 
on  the  rockets  in  a  CAMP ACT  container  irdergoing  a  longitudinal  impact. 
Results  ootained  from  this  mocei  should  not  be  considered  definitive,  but 
rather  should  indicate  whether  longitudinal  impacts  for  rocxets  inside 
CAMPACTs  should  be  more  or  less  severe  tr.an  impacts  on  unyieLding  surfaces. 


This  model  is  a  simple  extension  of  the  M55  FE  model  just  d 
the  nose,  a  nonlinear  elastic  element  has  Deer,  added  to  reoresen 
foam  insulation  or.  tr.e  CAM  PACT  door.  Also  at  the  nose,  a  lumped 
added  to  represent  tr.e  pallet  components  behind  the  modeled  M55 
weight  of  the  pa. let  containing  the  modeled  M55  rocket,  and  the 


t  the 


rocket 


and  roller  asserr.oiy  of  tr.e  comoact .  These  weights  are  added  at 
the  M55  rocket  '  rather  tr.an  tr.e  rear'i  so  that  their  loading  effe 
the  foam  only,  and  not  tr.e  rocket.  It  is  assumed  in  a  CAMP. ACT  i 
impact,  that  r.o  pallet  inertia  .cads  pass  through  the  rccwets.  o 
through  other  pal  let  mer.se  rs  ar.c  tr.en  into  the  foam.  This  ass  urn 
suoDorted  bv  tr.e  calcu.at  ions  of  Appendix  D. 
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2.2.3  Finite  Element  Character  1st  ics  ar.d  Material  Models 

The  M55  casing  is  fabricated  from  steel  ar.d  aluminum.  From  drawings 
provided  to  SwRI.  required  minimum  mechanical  properties  ar.d  or  the  onemical 
composition  of  the  casing  materials  were  found.  From  these  data,  the  ager.t 
cannister  (aluminum)  was  taker  to  oe  o06l  in  tr.e  76  temper  condition,  and  the 
rocket  casing  (stee.)  to  ce  AISI  ‘030,  water-quenched  ar.d  tempered  at  t00;F. 
For  aluminum  or  steel  finite  elements,  a  bilinear  el ast i c-o last i c  ~ateriai 
model  was  used.  The  tangent  moduli  for  the  steel  ar  :  the  n uni  nun  were 
calculated  cased  upon  ”9?  and  J"  strain  at  failure,  respectively.  As 
previously  discussed,  this  value  represents  one-half  of  the  .niaxiai  tension 
test  value  for  the  aluminum.  Pertinent  material  mod-: .  input  data  are 
summarized  below. 


Mechanical  Properties  of  AISI  1030 
Longitudinal  Impact  FE  Model 


Value 


Elastic  Modulus,  E 
Poisson's  Ratio 
Yield  Strength 
Tangent  Modulus,  E 


30.0  E  6 
0.270 
50.0  E  3 
1.39  E  5 
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Elastic  Modulus. 
Poisson's  Ratio 
Yield  Strength 
Tangent  Modulus. 


10.0  E  6 
0 . 300 
38.8  E  3 
88.5  E  3 


The  truss  element  representing  a  portion  of  the  thick  foam  layer  on  the 
CAMPACT  door  was  developed  from  test  data  [6].  The  length  of  this  element 
approximately  equalled  the  foam  layer  thickness  and  its  cross  sectional  area 
equalled  1  60th  of  the  door  area.  (In  the  CAMPACT,  60  rockets  impact  the  ocor 
during  a  longitudinal  impact).  The  stress-strain  characteristic  of  the  foam 
is  essentially  piece-wise  linear  and  consists  of  three  pieces.  The  modulus  of 
the  elastic  region  is  defined  as  "E"  and  tangent  moduli  of  p.astic  regions  are 
defined  as  "El"  and  "E2",  respectively. 


Table  3. 


Mechanical  Prooerties  of  CAMPACT  Foam 
Longitudinal  Impact  FE  Model 


Prooerf, 


Elastic  Modulus, 
Yield  Strength 
Tangent  Modulus, 
(at  k*  strain) 
Tangent  Modulus, 
(at  527.  strain) 


Value 

35.0  E  2 
14.0  E  1 
0.0 

99.2  E  1 


2 . 3  Lateral  Impact 

2.3.1  Analyses  Performed 

SwRI  calculated  the  dynamic  response  of  a  pallet  wh<=n  falling  from  uo 
feet  in  a  lateral  orientation  and  impacting  an  unyielding  surface.  Dynamic 
response  of  the  agent  cannister  due  to  impact  was  assumed  to  consist  of  two 
separate  phenomenon:  bending  of  the  rockets  and  squeezing  of  the  rockets  and 
tube  assemblies  between  the  lateral  wooden  supports.  Strains  and  stresses 
induced  by  these  phenomena  *ere  calculated  separately  and  superposed  to  obtain 
the  total  strain  or  stress  in  the  agent  cannister. 

A  two-d imens lonai  n lane  stress  FE  model  of  the  rocket,  launching  tube, 
and  lateral  support  and  a  irree-dimensionai  model  of  the  middle  row  of  M55 
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rocKets  in  the  pallet  assembly  were  created  to  compute  the  squeezing  and 
bending  results.  Figures  4  and  5  illustrate  the  prir.cica!  Features  contained 
in  these  finite  element  models. 

The  2-D  squeezing  analysis  (Figure  4)  was  a  calculation  of  the  quasi- 
static  response  of  the  tuoe,  agent  cannister  and  support  under  a  steadily 
increasing  lateral  shear  load.  From  tr.is  analysis,  shear  failure  load  of  tr.e 
rocket  tube  assembly  at  the  support  could  oe  assessed.  Figure  c  shows  tne 
fiberglass  maximum  principal  stress  versus  tne  total  shear  .cad.  as  we. i  as 
the  aluminum  cannister  strain  versus  total  sr.ear  load.  3v  linear  mg  tr.e 
plots,  one  can  see  that  when  the  fiberglass  reaches  a  stress  of  2  *  -:si  (its 
failure  stress),  strain  in  tne  aluminum  agent  cannister  is  near.',  2* .  Inis 
failure  stress  ar.d  cannister  strain  correspond  to  a  tota.  sr.ear  load  of 
approximately  25  kips.  Alt'nougn  the  2-D  analysis  *as  cor.t.rued  for  shear 
loads  greater  than  25  <ips,  these  results  inc.ude  the  stiffness  of  tne 
fiberglass  ard  are  r.ot  very  meaningful. 

It  is  postulated  that  at  tr.e  25  kips  shear  load  ieve.  the  fiberg.ass 
failed  completely.  The  strain  in  the  aluminum  at  this  load  .evei,  nearly  ll. 
is  presumed  to  step  cnange  to  a  value  greater  than  4J,  upon  fail. .re  of  tr.e 
fiberglass.  Hence,  a  simultaneous  failure  of  the  agent  carmister  and 
fiberglass  launching  tube  is  presumed  to  occur.  For  this  reason,  25  kips  is 
adopted  as  tne  cannister  failure  shear  load  at  the  supccrt. 

The  three-dimensional  bending  model  (Figure  5)  was  formulated  to  give  the 
maximum  sending  strain  and  the  maximum  shear  load  in  the  rocket  cannister.  It 
was  anticipated  that  squeezing  ar.d  bending  components  «ou_d  be  superposed  to 
estimate  the  total  strain  state  in  the  agent  cannister.  Sucerposition  would 
be  performed  at  times  when  either  the  agent  cannister  tending  strain  maxima  or 
the  cannister  squeezing*  strain  maxima  occurred.  As  .noicated  in  the  results 
section,  this  superposition  was  not  ultimately  required  because  the  squeezing 
strain  component  at  tr.e  laterai  support,  by  itself,  exceeded  the  failure 
cr  iter  ia . 

3oth  a  weight  check  ar.d  an  eigenvalue  cr.ecx  were  made  :n  tne  3-D  pallet 
model.  The  purpose  of  these  check  runs  was  to  assure  tne  investigators  that 
the  model  had  the  correct  mass  and  stiffness  properties.  The  weight  of  the 
middle  row  of  rockets  and  lateral  supports  was  calculated  by  nan:  to  ce  147 
lbs.  Weight  of  the  finite  element  model  was  (46  lbs.  Similar  agreement  was 
obtained  between  closed  form  ar.d  numerical  ca.cu iat ions  of  e . ger.va . ues .  A 
first  bending  mode  of  the  rocket/tube  assembly  was  caiculat-'d  to  oe  279  Hz.  A 
similar  mode  was  calculated  numerically  as  24 1  Hz.  Closed  form  eigenvalue 
calculations  are  contained  in  Appendix  3. 

2.3.2  Model  Development 

As  previously  mentioned,  F  ig  .ire  4  illustrates  the  feat  _r  :s  tr.e 
squeezing  model  used  in  the  lateral  impact  ara.ysis.  The  deve.oced  finite 
element  model  consisted  of  approximately  ' 90  2-1  p i me  stress  elements  naving 
a  thickness  approximately  equal  to  the  latera.  support  th.c-.r.ess. 

For  lateral  squeezing  of  the  agent  canristen  and  1  aunor.ir.g  tube,  only  one 
plane  of  symmetry  exists,  tne  plane  containing  the  lateral  centerline  of  the 
rocket  -  tube  assembly.  Thus,  tne  entire  ootf on  half  of  tr.e  rocket  and  t  aoe 
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assembly  are  represented  in  this  model.  .r.ey  are  modeled  with  beam  ar.d  slate 
elements  so  that,  with  minima*  modification,  the  lateral  squeezing  model  could 
also  oe  used  in  the  vertica.  crushing  analysis.  It  is  not  necessary  to 
represent  in  the  model  tne  entire  cotton  half  of  the  wooaen  support.  Under 
the  lateral  shear  loads  indicated  in  Figure  4,  only  the  right  portion  of  the 
support  is  .oacec  and  tr.us  on. .y  t.ois  region,  is  represented. 

Boundary  conditions  at  the  edges  of  the  lateral  support ,  ar.d  the  rocket 
casing  and  laun.cn Lr.g  tune  centero.ar.es  are  as  indicated  in  Figure  .  in 
addition,  nodes  ai.org  tne  agent  cannister-iauncning  tuoe.  arc  tube-support 
interfaces  are  ccuDled  m  tne  radial  direction.  7 he  purpose  of  these 
constraint  eo. at  ions  is  to  a.,  low  circumferential  slippage  between  the  M55 
rocket  and  its  .auncr.ing  tuoe,  as  occurs  in  reality.  Beam  and  o^ate  elements 
representing  the  .a.ncning  t.ce  and  agent  car.n.ister  are  also  coupled  via 
constraint  equations .  These  ouations  rectified  the  oeam  rotational  with  the 
plate  trar.s.at ior.a.  degrees  of  freedom. 

Figure  5  illustrates  tr.e  features  of  the  3-D  bending  model.  Only  tne 
middle  row  of  M55 ' s  in  the  pa. let  assembly  are  represented.  Because  there  are 
no  significant  .cad  patns  between  rows  of  rockets  in  tne  pallet,  interaction 
between  rows  under  lateral  impact  is  neglected.  The  FE  model  accounts  for 
interaction  oetween  rockets  comprising  the  middle  row  oy  connecting 
rocket  tuPe  assemb.ies  witn  truss  elements  representing  the  grain-direction 
stiffness  of  tr.e  .aterai  wooden  supports. 

To  obtain  correct  stresses  and  strains  at  launching  tube  and  agent 
cannister  cress  sections,  rocket 'tube  assemblies  were  modeled  by  coupled  sets 
of  oeam  elements,  or.e  set  of  elements  representing  tne  tubes,  one  set 
representing  the  rockets.  Constraint  equations  were  -.opiied  that  required 
transitions  (in  the  imDact  direction  only)  of  roexe:  and  tube  be  identical. 
These  coupled  beam  sets  are  represented  by  a  single  .me  in  Figure  5. 
Derivations  of  the  constraints  are  given  in  Appendix  3. 

Only  naif  the  longitudinal  length  of  the  pallet  is  represented  by  the  FE 
model.  Symmetry  has  teen  assumed  aoout  the  pallet  lateral  centerplane;  the 
effects  of  asymmetrical  response  modes  under  lateral  impact  are  assumed  to  be 
negligible  at  the  agent  cannister.  Boundary  conditions  at  the  rockets' 
midpoints  are  as  sr.own  in  Figure  5.  Launching  tunes  had  identical  boundary 
conditions  at  corresponding  nodes.  Elements  representing  the  rockets  had 
annular  cross  sections  and  densities  which  accounted  for  casing  and  internal 
weigr.ts  such  as  the  chemical  agent  and  solid  roexet  fuel. 

2.3.3  Finite  Element  Character  ist  ics  and  Material  Models 

The  squeezing  and  bending  materials  developed  for  tr.e  lateral  impact 
problem  were  identical  except  tr.at  or.e  additional  material  model  existed  in 
the  3-D  model  that  did  not  exist  m  the  2-3  analysis  (AISI  "030  stee.).  As  m 
the  case  of  the  longitudinal  impact  problem,  aluminum  and  steel  rocket  casings 
were  represented  as  elast ic-plastio  materials  with  tangent  moduli  calculated 
based  upon  failure  strains.  Again,  for  the  aluminum,  the  value  used 
represented  or.e -half  tf  the  mi  axial  tension  test  data.  Pertinent  material 
model  input  data  are  summarized  sol 
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Mechanical  Properties  of  A. SI  1030  Steel 
lateral  Impact  FE  Model 


Elastic  Modulus. 
Poisson's  Fatio 
field  Strength 
Tangent  Modulus, 


30.0  E  6 
0.270 
90.0  E  3 
1.39  E  5 


Table  5. 


Mechanical  Properties  of  6061-T6  Aluminum 
Lateral  Impact  FE  Model 


Proper  r> 


Elastic  Modules, 
Poisson's  Fatio 
field  Strength 
Tangent  Modulus, 


10.0  E  6  psi 
0.300 

38.3  E  3  psi 


Fiberglass  material  data  were  developed  from  the  tube  specif ication  given 
in  the  pallet  assembly  drawings  and  from  Owers-Corning  Fiberglass  Company  data 
found  in  reference  7.  An  average  glass  to  resin  ratio  of  355  was  assumed.  A 
sufficient  amount  of  this  glass  was  assumed  to  be  in  mat  form  to  resuit  in 
isotropy . 

Fiberglass  is  a  brittle  elastic  material,  i.e.,  there  is  1  linear 
relation  between  stress  and  strain  up  until  failure;  however,  material  models 
with  ultimate  strength  cut-offs  do  not  exist  in  the  3wPl  version  of  AETNA. 
Hence,  for  the  purposes  of  computation,  an  elastic-perfect ly  clastic  material 
model  was  used  for  the  fiberglass  launching  tubes.  Following  calculation  of  a 
solution,  launching  tube  elements  were  manually  cnecked  for  "elasticity" 
(failure).  For  the  40  foot  drop  height,  "piast ic i ty "  in  the  fiberglass  was 
insignificant,  occuring  only  over  a  small  portion  of  the  tube  length  and  never 
occurring  through  the  complete  wall  thickness. 


Mechanical  properties  used  for  the  launching  tubes  are  summarized  below. 


7aoie  6. 


Mechanical  Properties  of  Fiberglass  Launching  Tunes 
Lateral  Impact  FE  Model 
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Axial  or  grain-direction  stiffness  of  the  wooden  rocket  supports  was 
estimated  from  an  effective  cross  sectional  area  for  the  support .  Because 
supports  have  a  "scalloped"  shape,  their  actual  area  varies  along  the  support 
length.  To  simplify  the  modeling  of  these  members,  an  effective  team  width 
was  computed  based  on  the  oeam  longitudinal  area  and  lengtn.  (See  Appendix 
B.)  Because  most  of  the  stresses  during  lateral  impact  were  expected  to  be 
axial,  the  beam  element  in  this  model  was  given  isotropic  material  properties 
identical  to  the  orthotropic  grain  properties.  Yield  strength  was  taken  at 
12*  moisture  content,  conforming  to  seasoned  wood  [8],  [9]. 

Pertinent  character isc ics  of  these  elements  are  shown  in  the  table  below. 


Table  7. 

Mechanical  Properties  Wooden  Lateral.  Supports 
Lateral  Imoact  Modei 


Property 

Elastic  Modulus. 
Poisson's  Ratio 
Yield  Strength 
Tangent  Modulus, 


Value 

1.353  E  6 
0.239 


5700.0 

0.0 


03  1 


03  1 
03  1 


2 . 4  Vertical.  Impact 

2.4.1  Analyses  Performed 

SwRI  computed  the  resoor.se  of  the  pallet  assembly  when  failing  40  feet 
onto  a  unyielding  surface,  in  its  normal,  upright  orientation.  Dynamic 
response  of  the  agent  cannister  to  impact  was  assumed  to  consist  of  two 
independent  phenomena:  bending  of  the  rockets  between  lateral  supports,  and 
diametral  crushing  of  the  iauncr.ing  tune  and  rocket  cross  sections  cy  the 
supports.  Strains  and  stresses  induced  by  these  phenomena  were  computed 
independently  and  superposed  to  obtain  the  total  strain  or  stress  in  the  agent 
earn  ister . 
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■.-?  invest  .ga.crs  created  :  mice  element  mea*  izat  ior.3  to  compute  ter.dir. 
r\.sr.ir.i  results .  A  -c  i.mer.sionai ,  plane  stress  model  of  the  rocket . 
ting  :.:e,  and  support  interface,  and  a  three-c imens  tonal  model  of  one 
-jr  of  tr.e  pallet  assert*-;  were  -.sea  to  obtain  crushing  and  bending 
ses.  Figures  "  ana  3  illustrate  the  mam  features  of  these  models,  and 
::x  det  a.  .s  tr.e  effective  density,  concentrated  mass,  and  team  propert 
.  it.cns  ”.tde  in  creating  the  FE  models. 

7ne  :r. sting  analysis  was  a  calculation  of  the  quasi-static  response  of 
,te.  agent  cannister,  and  wooden  support  uncer  a  monotonicaily  increas.r. 
.ng  load.  From  these  computations,  a  crushing  force-cam ister  strain 
ctenstic  was  developed.  Initial  conditions  for  the  dynamic  response 
s.s  tf  tr.e  three-dimensional  motel  corresponded  to  the  40  foot  drop 
results  from  this  model  were  strain-time  histories  in  the  roe-cet 
c  sue  to  tending,  and  the  crushing  force-time  history  at  the  rocket- 
i-  supoort  connections. 

Toe  crusting  force-cam  ister  strain  characteristic  was  cotained  first, 
tr.ese  data,  a  nonlinear  elastic  truss  element  ("tr.e  crushing  spring",'  wa 
:tec  wrier,  reflected  the  computed  roctcet/ tube  and  lateral  support 
mg  st.ffress.  These  truss  elements  became  part  of  the  3-D  bending 
Hence,  tr.e  dynamic  response  computed  using  the  3-D  model  reflected 
tending  and  crushing  energy  aosorption. 

Both  a  weignt  checx  arc  an  eigenvalue  check  were  made  on  the  3-D  pallet 
tefere  computing  dynamic  response.  The  purpose  of  these  calculations 
c  assure  tr.e  investigators  that  the  model  had  the  correct  mass  and 
ness  properties.  Tr.e  weight  tf  the  quarter  pallet  modeled  was  estimated 
lit  Ids.  A  we.ght  tf  126.4  Ids.  was  calculated  from  the  FE  model.  Ver 
agreement  was  a. so  obtained  between  closed  form  and  numerical  solutions 
iger.values .  A  first  censing  mode  cf  the  rocket  tube  assembly  was 
-atec  to  ce  y.»  Hz.  A  similar  rocket  -'tube  bending  mode  was  calculated 
ica.. y  as  293  Hz.  Eigenvalue  calculations  are  shown  in  Appendix  C. 

1.-.2  Model  Pevelcomert 

To  obtain  the  crusnir.g  force-cannister  strain  characteristic ,  a  finite 
rt  mcdel  tf  the  rocket,  tuoe  and  support  interface  was  created.  Because 
r.g  arc  geometrical  symmetry  exists  about  the  vertical  and  horizontal 
rp.ar.es  of  this  connection,  only  one-quarter  of  the  tube  and  support 
n  need  be  modeled.  Figure  7  shows  the  principle  features  of  the  finite 
r.t  mode*  developed.  This  idealization  consisted  of  about  ’00  2-D  plate 
r.ts  having  2  DOF  at  each  node.  The  thickness  of  these  elements  equalled 
nickness  of  the  wooden  lateral  support. 

Fo.r.dary  conditions  at  the  edges  of  the  lateral  support  member,  and  the 
t  casing,  and  launching  tube  centerplanes  are  indicated  in  Figure  7.  In 
.or,  nodes  along  the  agent  car.n  ister- launch  mg  tube  and  tube-support 
faces  are  coupled  ir.  the  radial  direction.  These  constraints  allow 
mferent.al  s.ippage  between  the  M55  and  its  launching  tuoe  as  actually 


F : gure  3  shows  the  features  of  the  bending  model.  One-quarter  of  the 
t  containing  three  co**r.rs  of  rockets  and  tubes  are  represented.  The 
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off-center  columns  of  rockers  are  idea* ized  is  lumped  masses  on  the  latera* 
support  teams,  while  the  center  column  cf  rccxets  are  mere  explicitly  mcae*ed 
using  beam  elements.  "Crusr.ir.g  SDring"  elements  are  indicated  in  me  figure. 


In  order  to  account  for  possicie  contact  oetween  the  impact  surface  and 
the  lateral  support,  a  truss  element  n.as  oeen  installed  at  the  support  center 
point  (see  Figure  3).  This  element  had  little  axial  stiffness  for 
displacements  less  than  the  support-surface  gap;  for  greater  displacements, 
the  element  stiffness  increased  sharply. 

Again,  only  half  the  longitudinal  length  of  the  pallet  is  represented  by 
the  model.  Symmetry  ras  been  assumed  aoout  the  pallet  lateral  centerpiar.e, 
and  the  effect  of  asymmetrical  modes  under  vertical  impact  are  assumed  t:  be 
small  at  the  agent  cannister. 

2.4.3  Finite  Element  Character  1st ics  and  Material  Models 

The  3-D  model  developed  for  the  pallet  vertical  Impact  proclem  was  me 
most  complex  of  the  pallet  models  constructed.  This  idealization  of  or.e- 
quarter  of  the  pallet  consisted  of  over  '00  beam  and  truss  elements  classified 
into  six  different  groups,  each  group  having  particular  dimensional  or 
constitutive  features.  In  the  following  paragraphs,  each  group  is  discussed 
in  turn. 


Lateral  wooden  supports  in  the  pallet  were  idealized  as  beams  of  uni 
cross  section  ana  nomoger.ous.  Isotropic  material.  Because  lateral  suppor 
have  a  "scalloped"  shape  so  as  to  conform  to  the  laurehing  tubes,  their  ; 
cross  sectior.ai  area  varies  along  the  support  length.  To  simplify  the 
modeling  of  these  mencers  while  accounting  for  the  regions  of  increased  : 
sectional  area  between  "scallops",  an  effective  beam  width  was  computed  : 
upon  the  actual  beam  Longitudinal  area  and  length  (see  Appendix  C).  Bee; 
most  of  the  stress  was  expected  to  be  caused  by  support  bending  along  tr.t 
grain  direction,  this  beam  was  given  isotropic  material  properties 
corresponding  to  the  grain  direction  orthotropic  properties.  Because  po; 
yield  cnaracter istics  of  woods  are  r.ot  well  understood,  it  was  assumed  tr 
the  supports  had  eiast ic-per feet ly  plastic  oehavior.  Yield  strergtn  for 
support  was  taken  for  wood  at  12^  moisture  content  (seasoned  wood).  Pol; 
ratio  was  taken  as  the  average  of  the  radial-grain  and  tangential-gram 
ratios . 

Pertinent  characteristics  of  these  elements  are  summarized  in  the  r.t 
table . 
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Table  3. 

Mechanical  Properties  Wooden  Lateral  Supports 
Vertical  Impact  Model 


Property 

Elastic  Modulus, 
Poisson's  Ratio 
Yield  Strength 
Tangent  Modulus, 


Value 

1.353  E  6 
0.239 
5700 . 0 
0.0 


.he  wooden  base  stringer  .and  ruccir.g  strics  a, erg  tr.e  ret  tom  o:  tr.e 
pallet  were  assumed  to  sustain  only  cross  gram  corner  ess  ive  .ties  during  the 
impacts.  Thus,  like  the  lateral  supports,  these  ta; let  ccrpcrer.ts  were 
idealized  as  isotropic,  e Last  ic-cer feet  1 plastic  mater :a.s  -v.tr.  to:  .h  ard 
mechanical  properties  correscorc.r.g  to  tr.e  cross-zrair  wood  crthotroc ic  axis 


Parameters  input  :or  tr.ese 


:ts  are  sr.cwr. 


.  ao.e  i. 


Mectar ica„  Properties  r  or  ease  ctrmger  are  rubcir.g 
Vert  ioa .  ”0  MO20I 


rrcp0r tv  : 

Elastic  Modulus,  Ep  O.Ot 
Yield  Strength  ‘  510.0 
Taneer.t  Modulus,  E~.  0.0 


As  in  the  case  of  the  Longitudinal  and  .aterai.  impact  problems,  rocket 
casing  materials  were  considered  to  be  isotropic,  elast ic-piast ic  materials 
with  the  following  properties. 


Taole  10. 

Mechanical  Properties  of  AISI  10] 
Vertical  Imoact  EE  Model 


Property 


Elastic  Modulus, 
Poisson's  Ratio 
Yield  Strer.gtr. 
Tar.gert  Modulus, 


30.0  E  5 
0.2  70 
90.0  E  3 
1.39  7  : 


Mechanical  Properties  of  6o6'-T6  Aluminum 
Vertical  Imoact  EE  Mocel 


Procer t\ 


Elastic  Modulus,  E 
Poisson’s  Ratio 
Yield  Strength 
Tangent  Modulus,  E‘ 


‘0.0  E  6 
0 . 300 
33.3  E  3 
53.0  E  3 
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found  in  reference  o.  An  average  glass  to  ream  ratio  or'  35  • 
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efficient  amount  cf  this  glass  .»i.s 
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isotropy.  Fiberglass  is  a  crioo.e  elastic  mater  :ai .  i.e.,  tr.ere  is  j  mr.-ir 
relation  between  stress  and  strain  „p  ur.t : .  r'i.i.ra.  Materia,  me  tel  a  „  it", 
ultimate  strength  out-offs  do  ret  exist  m  me  2.mi  /ersior.  of  I  i '  1 A .  Feme, 
for  computatiors.  an  ei as: ic-cer fee:  .v  p.sstic  -rater  ia.  -.ode.  as  .see  for  th 
f ioerglass  launching  tunes,  Fol. owing  cal  cu.  at :  or.  of  i  so.ut.cn.  1  i.-cr.  mg 
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drop  r.eigr.t,  no  "plasticity"  occurred. 
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Mechanical  Frooerties  of  F i : 


.  x  .  1*0*55 


Elastic  Modulus, 
Poisson’s  Patio 
Yield  Strength 
Tangent  Modulus, 


’  .  0  E 
0.  ’  ■ 

2  .  4  0  E 
0.0 


The  non.  mean  elastic  truss  element  represent  in  •  tr.e  diame 
stiffness  of  tr.e  cannister  cross  section  was  develot-.:  from  tr.e 
cata  computed  .str.g  the  2-D  mocel.  Displacement  at  i ,  indicate 
and  the  total  applied  force  /.ere  tabulated  from  the  2-D  analyse 
uere  converted  into  an  eq._va.ent  stress-strain  character  ist  ic 
element  length  and  cross  seccicrai  area  of 

3.03  inches  and  1  square  inch,  respect  ivei y .  Figure  3  sr.ows  th> 
cispiacenert  data  cbta.ned  from  tne  2-D  results .  The  data  com: 


trai  crushmu 
cr-sr. irg  f:r 


.  r.ese  1-: 
as  sum  ing 


inch,  F  =  7200  lbs .  corresponds  to  failure  at  the  agent  cannist 
strain  through  the  wa.l  of  tr.e  cannister.  Appendix  2  contains 
me  displacement -maximum  cannister  strain  results. 
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Lateral  supports  in 
stringers,  also  fabricate 
impacts,  these  lor.gitudin 
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Mechanical  Properties  c:  uong itud ir.ai 
Vertical  Imcact  FE  Model 


Elastic  Modulus, 
Yield  Strengtn 
Tangent  Modulus, 


Value 

0.064  E  6 
510.0 
0.0 


Under  static  loads,  such  as  gravity,  the  pallet  assembly  is  sucport 
entirely  by  tr.e  large  base  stringer  ar.d  rubbing  strip.  Under  impact  lea 
however,  inertia  loadings  may  ce  sur'f icientiy  high  so  that  lateral  succor' 
teams  undergoing  tending  ai.so  contact  the  impact  surface.  To  represent  tr 
change  ir.  boundary  conditions  ordinarily  calls  for  a  nonlinear  element 
specifically  adapted  to  contact  prooiems.  in  SwRI's  version  of  A31NA,  sue 
element  ices  not  exist.  To  address  this  difficulty,  tr.e  investigators  ere 
a  nonlinear  elastic  truss  element  with  increasing  stiffness.  This  element 
shown  in  Figure  3,  was  installed  at  the  midpoint  of  the  lateral  support, 
support  displacements  less  than  cr  equal  to  the  distance  from  the  suopert 
the  impact  surface,  this  element's  stiffness  was  very  low.  For  large 
displacements,  the  element’s  stiffness  increased  rapidly,  representing  sui 
contact  with  the  impact  surface. 
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As  mentioned  previously .  S*?l  considered  three  ixpact  scer.ar :c 
caused  'ey  drops  fren  -0  ar.d  30  feet  onto  rigid,  unyielding  surfaces 
ixpact  scenario  tr.at  considered  tr.e  additional  mass  ar.d  energy  aese 
materials  .  foam)  of  tr.e  1  AM?  ACT.  In  addition,  the  possibility  of  f 
the  agent  cam  ister  oy  ouc.e.mg  was  also  examined. 


ar.d  an 


r  inures 


surr.ari.de  tr.e  results  of  the  longitudinal  anyiel: 


surface  impact  analyses,  figure  *0  shows  the  locations  a.or.g  tr.e  agent 
cann ister  wnere  axial  strain  data  were  taou^ated.  Figure  *  i  contains  three 
plots  of  the  strain-time  r.i story  for  the  -*0  foot  impact.  7 he  plots  shew 
strains  at  the  forward,  middle,  and  aft  er.d  of  the  agent  cann  ister . 

Similarly,  Figure  '2  shows  o-ots  of  the  strain-time  history  for  the  SO  foot 
impact  case.  It  consists  of  plots  of  the  strain  response  at  6  locations  a.  or  2 
the  agent  cannister. 

For  the  40  foot  drop,  the  strains  in  Element  5  (topmost  plot  of  Figure 
II)  exceeded  the  agent  cannister  failure  criterion  by  a  good  margin,  reaching 
8%  at  .0035  seconds  after  impact.  These  high  strains  occurred  in  the  forward, 
tapered  region  of  the  cannister.  At  points  aft  of  the  tapered  region  where 
cross-sectional  areas  are  larger,  lower  strains  were  recorded.  At  Elements  11 
and  15,  maximum  strains  computed  were  0. 35?  and  0.1%,  respectively. 

Close  examination  of  these  plots,  particularly  tr.e  topmost  plot  of  Figure 
“1,  reveals  flat  regions  in  the  strain-time  history.  These  flat  response 
regions  are  generally  preceded  by  a  rapid  increase  ir.  strain  and  only  occur 
after  material  strains  have  exceeded  0.4%  (yield  strain). 

The  rapid  rise  in  strain  before  the  flat  regions  signifies  that  the 
applied  load  is  rising  ar.d  that  material  yielding  is  ocurring;  the  strain 
increases  rapicly  because  the  plastic  stress-strain  modulus  is  very  low 
compared  to  tne  elastic  modulus.  After  the  load  begins  to  decrease,  however; 
the  material  response  can  again  be  elastic.  If  the  applied  load  increases  for 
a  second  time,  the  change  in  strain  will  be  relatively  small  until  the  new 
yield  point  is  exceeded.  These  flat  spots  in  the  strain  resDonse  represent 
points  in  the  strain  time  history  where  this  elastic  behavior,  below  a  r.ewiy 
defined  yield  point,  is  occurring. 

An  equivalent  crush  force  for  impact  conditions  was  computed  for  this 
scenario  by  equating  the  pallet's  kinetic  energy  at  impact  to  the  work 
necessary  to  deform  the  structure.  For  longitudinal  impact,  this  equivalent 
crush  force  was  546.800  ibs.  (see  Appendix  E).  Note  that  failure  of  the 
rocket  occurred  in  this  impact  scenario. 

For  the  30  foot  impact  analysis,  strains  at  6  points  along  the  agent 
cannister  were  plotted  (Figure  ’2).  Again,  as  in  the  40  foot  imDact  case, 
strains  were  greatest  in  tne  forward,  tapered  regions  of  the  cannister  and 
decreased  at  aft  sections.  Strains  at  ail  sections  were  consistently  lower 
than  the  40  foot  case,  as  expected.  Maximum  strain  at  the  foremost  cannister 
section  was  less  than  1.5%,  indicating  that  yielding  of  the  casing  had 
occurred,  but  it  was  not  sufficient  to  cause  failure  (yie^d  in  aluminum  begins 
at  approximately  0.4%  strain).  Maximum  strains  at  aft  sections  were 
approximately  0.31%,  0.30%,  0.29%,  and  .235%,  respectively.  Strains  at  the 
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(eont'd)  Axial  Strains  in  M55  Rocket  Undergoing 
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Mote  that  strains  caused  ov  30  ar.d  -0  foot  impacts  are  r.ot 
preport  ionai .  Crops  from  30  feet  cause  approximately  :r.iy  ar.e-haif  the  stra: 
in  the  forward  canmster  caused  tv  tr.e  40  foot  impact.  Impact  responses  of 
the  earn  ister  are  rot  proccrt  ioral  cecause  the  canr.ister  respor.se  Is  non¬ 
linear  .  Agent  canmster  elements  5  through  3,  for  example,  deformed 
plastically  curing  tr.e  -0  foot  imoaet,  while  during  tr.e  30  foot  impact  or.iv 
car.nister  element  5  deformed  plastically  -  the  greater  part  of  the  agent 
earn ister  did  ret  vie.d. 


Figures  13  ar.d  summarize  the  results  for  the  foam  impact  model.  As 
previously  discussed,  tr.e  model  attempts  to  assess,  aloeit  in  a  qualitative 


way,  interactions  cetweer.  the  oal lets  and  the  CAMPAC7  structure 
longitudinal  impact. 


it  appears  that  the  irner-iiner  tf  the  CAMPACT  has  the  potential  to 
significantly  ameliorate  tr.e  response  of  the  pallets.  Figure  13  shows  that, 
for  longitudinal  impacts  onto  tr.ick  foam,  the  M55  rocket  response  is  primarily 
in  one  mode  (compared  to  tr.e  multiple  mode  response  indicated  in  Figures  1‘ 
and  '2  for  unyielding  surface  impacts)  and  that  the  maximum  strain  response  is 
substantially  reduced.  Maximum  computed  dynamic  strain  for  the  agent 
car.nister  in  this  analysis  occurs  in  the  forward  tapered  section  of  the 
cannister  and  was  aoout  .0415  as  compared  to  85  strain  for  the  unyielding 
surface  impact  case.  (Because  structural  damping  was  not  included  in  these 
analyses,  strain  response  amplitudes  do  not  decrease  over  time). 


Besides  failure  oy  ovei ,^ad,  another  possible  failure  mode  for  the  agent 
car.nister  is  buckling  oetween  the  lateral  supports.  Four  bucxiing  moces  were 
examined  in  order  to  assess  whether  buckling  was  a  ci it  leal  failure  mode  for 
this  type  of  impact.  In  the  first  case,  the  cannister  was  considered  as  a 
Euler  column  and  the  critical  axial  load  was  computed.  Because  the  agent 
cannister  is  a  thin  walled  cylinder,  in  the  second  case  lobar  or 
circumferential  ruckling  was  assessed.  In  this  analysis,  the  agent  cannister 
was  considered  to  have  uniform  cross-sectional  area,  i.e.,  the  forward  tapered 
section  was  ignored.  In  tr.e  third  case,  lobar  buckling  of  the  taper  section 
was  assessed.  Fere  ar.d  aft  cross  sections  of  the  taper  section  were  presumed 
to  remain  circular  in  this  closed  form  solution.  In  the  last  case,  the 
internal  burster  casing  critical  buckling  load  was  computed,  again  assuming 
that  the  burster  acted  as  a  Euler  column.  Results  of  the  bucxling  analyses 
are  summarized  in  Table  "4.  in  all  cases,  the  probable  buckling  stress 
exceeds  the  agent  cannister  ultimate  stress  (42.0  ksi).  Hence,  for  these 
rocket  support  conditions,  failure  by  overload  should  occur  before  failure  by 
buckling. 


'.g  i  cud 


esults  indicate  that  a  high  probability  of  agent  leakage  exists 
inai,  unyielding  surface  impacts  from  heights  equal  to  or  greater 
Pallet  assemblies  impacting  unyielding  surfaces  from  lower 
e  a  much  lower  probability  of  catastrophic  failure  of  the  agent 
The  results  also  appear  to  indicate  that  the  CAMPACT  assembly 
liorates  the  effects  of  longitudinal  impact  significantly, 
hout  a  mere  ccm.pi.ete  cal  let  CAMPACT  interaction  analysis,  it  is 
to  say  what  the  safe  net  safe  drop  height  might  oe  for  M55 
ced  in  tr.e  CAMPACT. 
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Figures  *5--7  shew  tr.e  results  computed  far  tr.e  .aterai  impact 
scenario.  Figure  '5  illustrates  tr.e  locations  along  tr.e  M55  roexet  where 
strain  data  from  the  3-D  analysis  were  obtained.  These  data  aDpear  as  tr.e 
slots  of  Figure  '6.  E.erer.ts  58  tnrougn  -*3  represent  the  agent  car.nister 
region  of  the  rocket,  with  element  -*3  representing  its  forward  tapered 
section.  Figure  '7  snows  tr.e  roexet  tube  assemoly  snear  force-time  motor 
for  tr.e  forward  ar.c  aft  support  junctions. 


rate  tr.ai 


The  strain  respor.se-t  ime  mstories  oiotted  in  Figure  ’t  are  primarily  t  hi 
result  of  rocket  bending  net  ween  the  supports.  The  plots  indicate  tr.at  the 
impact  was  not  sufficient  to  cause  failure  of  the  cannister  at  locations  a*,  ay 
from  the  support  junctions.  The  maximum  bending  strain  observed  occurred  in 
Element  39.  aocut  2  mii.isecords  after  impact.  This  strain  exceeded  21,  or 
about  one -ha if  the  failure  value.  At  ail  other  elements  alo-g  the  cannister 
region,  oendir.c  strains  were  at  all  times  less  than  1 1. 


Squeezing  of  the  rocket. tube  assembly  between  the  latera.  supports  is 
quite  severe  in  this  loading  scenario,  however.  Figure  17  depicts  tr.e 
rocket. 'tuoe  shear  force  at  the  support  junctions  indicated  in  Figure  15.  The 
maximum  shear  force  occurs  around  2  milliseconds  after  initial  impact  at  the 
aft  lateral  support  junction.  "This  shear  force  momentarily  exceeds  50  xips. 
From  the  quasi-static  squeezing  analysis,  it  was  determined  that  a  25  kip 
shear  force  was  sufficient  to  cause  cannister  failure,  i.e.,  strains  exceeding 
of  41  (see  Sect  ion  2 . 3.  '  • .  Hence,  it  appears  that  failure  of  the  agent 
cannister  is  likely  to  occur  c  arm*  lateral  impacts . 


Again,  ty  equating  the  work  r.ecessarv  to  deform  tr.e  structure  to  the 
kinetic  energy  of  tr.e  pa._et,  an  equivalent  crush  force  was  computed  for  trie 
impact  scenario.  This  equivalent  force  equalled  627.950  lbs.  (see  Appendix 
E).  Note  that  m  this  scenario,  the  missile  failed  by  localized  crushing  cf 
the  azen.t  cannister  at  the  lateral  succor:. 


3.3  Vert  .cal  Impacts 


Figures  '5-20  shew  the  computed  results  for  tr.e  pallet  assembly 
undergoing  a  -0  foot  drop  onto  a  unyielding  surface.  Figure  '3  indicates  the 
approximate  locations  of  the  elements  for  which  results  „ere  plotted. 

Elements  21  through  26  model  tr.e  agent  car.nister,  with  Element  2 1  represent  ir.g 
the  forward,  tapered  region.  Element  27  has  the  properties  of  tr.e  rocket 
motor  casing.  Figure  '9  contains  the  plots  of  the  strain-time  histories 
obtained  from  the  3-D  finite  element  model,  and  Figure  20  summarizes  the 
diametral  crusr, ing  displacements. 


itrairs  ar 


:uite  .ow 


vne  mode  cf  response  of  tr.e  M55  rockets  during  a  vertical  impact  is  in 
tending ,  ar.c  it  is  apparent  from  Figure  '9  that  the  magnitude  of  tr.e  bending 
strains  are  quite  low.  These  low  strains  indicate  that  the  impact  energy  was 
not  ce  ing  absorbed  significantly  cv  rocket  bending.  (It  „.as  for  this  reason 
that  the  time  integration  analysis  were  stopped  at  4.5  milliseconds),  Fesu.ts 
from  the  3-D  analysis  show  that  significant  strains  and  yielding  occurred 
during  compression  of  the  base  stringer  rubbing  strip  ar.d  bending  cf  the 
.  iter  a  1  wooden  supports.  ?1  tsticity  «as  reserved  in  both  the  rubbing  strip 
and  lateral  supports.  The  maximum  c  i  sc  1  acemer.t  of  the  wooden  lateral  supports 
~  as  less  than  two  incr.es.  Th  _  3 .  it  joes  not  appear  that  catastrophic  fa.!  -re 
;  f  "e  supports  wou.d  occur  after  Jef.ect ions  of  this  magnitude.  Further, 
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Figure  19.  Strains  in  M55  Rocket  Undergoing 
Impact  (from  40  feet) 
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’.13  uer.eetion  ls  not  surriciert  to  oat 
-.tact  surface. 


Crushing  at  the  rocket  tube-supper 
cause  failure  in  the  agent  car.nister. 

:s  not  exoected  for  t.nis  leasing  sce.nar 


motion  also 


is  not  expected  for  t.nis  loading  scenario.  Figure  20 
displacements  at  the  forward  and  aft  support  junctions 
rocket.  (Fiat  regions  in  tr.ese  plots  were  caused  tv  r: 
plotting  of  the  results.  Actual  displacement  variat  :o: 
but  were  quite  small  ' .  Maximum  displacement  occurs  at 
equals  0.02  incr.es,  approximately .  From  the  quasi-sta: 
it  was  calculated  that  diametral,  crusr.ing  displacement: 
were  necessary  to  induce  car.nister  failure,  i.e.,  cam: 
than  •>%.  The  suoerpos it ior.  of  bending  ana  crushing  st: 
strain  much  less  tr.an  - u .  Crerefore,  it  is  considered 
is  uni  ike. y  to  oe  caused  by  this  impact. 


3.h  Cone  1  .d ir.g  Remarks 

Overall,  it  should  be  noted  that  tr.e  analysis  methodology  is  sor.ewh 
conservative,  particularly  because  of  the  assumption  that  impact  surface 
perfectly  rigid.  Actual  surfaces,  and  the  CAMPACT  structure  also,  will 
impact  energy,  thereby  reducing  car.nister  strains.  Conversely,  it  must 
forgotten  that  the  rocket  structure  was  idealized  in  these  analyses:  ag 
effects,  corrosion,  casting,  or  welding  flaws  are  not  assessed  or  accoun 
for.  3ecause  these  latter  effects  will  tend  to  induce  leakage,  .at  lower 
expected  drop  heights,  results  presented  here  must  oe  :sed  carefully,  as 
in  the  light  of  past  experience,  and  compared  with  actual  test  Java  w her 
ooss ib  le . 
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